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Universal spin-triplet superconducting correlations of Majorana fermions 
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We establish that Majorana fermions on the boundary of topological superconductors have only 
spin-triplet superconducting correlations independent of whether the bulk superconducting gap is 
spin singlet or triplet. This is universal for time-reversal broken (respected) topological superconduc¬ 
tors (TSCs) with an odd number of (pairs of) Majorana fermions on the boundary. Consequently, 
resonant Andreev reflection induced by Majorana fermions only occurs in spin-triplet channels and 
always injects spin-triplet Cooper pairs into the leads. This spin-triplet condensate results in a 
spin-orbit coupling (SOC) controlled oscillatory critical current without a 0-7T transition in the 
TSC/SOC-semiconductor/TSC Josephson junction. The observation of this unique current-phase 
relation can serve as the definitive signal for Majorana fermions. Our study shows a technique for 
manipulating Majorana fermions based on their spin-triplet superconducting correlations. 

PACS numbers: 74.45.+C, 75.70.Tj, 85.25.Cp 


Introduction- Due to recent theoretical proposals 
of realizing topological superconductors (TSCs) with 
an s-wave superconducting proximity effect in spin- 
orbit coupling (SOC) systems 1-13 , finding and ma¬ 
nipulating Majorana fermions 14-16 have become possi¬ 
ble in superconductor-semiconductor devices 17-23 . Ma¬ 
jorana fermions (MFs) in this context are symme¬ 
try protected zero-energy subgap states of TSCs with 
equal weight electron-hole superposition. This leads 
to some unique transport properties, such as quan¬ 
tized zero-bias conductance 24-26 and fractional Joseph¬ 
son effect 5,10,15,27-29 , and exotic superconducting con¬ 
densates such as an odd frequency 8-wave and spin- 
triplet pairing in a time-reversal broken 77-wave- 
superconductor/normal-metal junction 30 . However, so 
far, most studies ignore the spin states of MFs. On the 
other hand, controlling or manipulating spin states is a 
conceptional revolution in fundamental science as well 
as emerging technologies 31 and is triggering important 
fundamental discoveries of spin-dependent transport phe¬ 
nomena such as giant magnetoresistance 32,33 and spin 
Hall effect 34,35 . Therefore, it is natural to ask whether 
we can manipulate MFs by their spin properties. 

In this paper, we show that MFs, zero-energy bound 
states of TSCs, have only spin-triplet superconducting 
correlations which are model independent and guaran¬ 
teed by particle-hole symmetry for time-reversal bro¬ 
ken TSCs (D class 36 ) and particle-hole as well as time- 
reversal symmetries for time-reversal respected TSCs 
(Dill class 36 ). This leads to our experimental predic¬ 
tion that no matter whether the bulk superconducting 
gap is spin singlet 7,25 or triplet 26,37,38 , the MF-induced 
resonant Andreev reflection (AR) only injects spin-triplet 
Cooper pairs into the normal lead. We perform a sym¬ 
metry analysis as well as numerical calculations of the 
reflection matrix at the interface of either time-reversal 
broken or respected TSC/normal-metal (NM) junctions. 
We further propose to detect spin-triplet Cooper pairs in 
a TSC/SOC-semiconductor/TSC experimental setup by 


observing a SOC controlled oscillatory critical Josephson 
current. It has been recently shown that SOC can rotate 
the d vector of the spin-triplet Cooper pairs and leave 
spin-singlet Cooper pairs unaffected in the proximity re¬ 
gion. Therefore, the observation of our predicted SOC 
controlled critical current can serve as direct confirma¬ 
tion for spin-triplet Cooper pairs and the MF transport 
signal. 

Spin-triplet superconducting correlations of MFs In 
any dimension superconductor, the MF, a self- 
Hermitian quasiparticle, can be defined as 7 i(x) = 

^7/7,7/7, ^{, — where x is generally a vec¬ 

tor and ^i(cc) is the ith zero-energy quasiparticle wave 
function. Electron and hole correlations at the Fermi 
surface show up in the block off-diagonal part of the 
equal space zero energy spectral function, which is a 
Hermitian matrix and has the form A(E = 0, as) = 
JA |^(a?))(J/i(a?)| 39 . In the four-component spinor ba¬ 
sis ^7/7, 7/45 Vfp — 5 the block off-diagonal part of the 
spectral function (BOSF) can be written as 


A oE (E,x) = 


0 

[(docro + d • cr )] f 


(d 0 cr 0 Pd a) 


where <ro and a are the 2x2 identity and Pauli matri¬ 
ces in electron-hole and spin space respectively, and do 
and d are the spin-singlet and spin-triplet pairing am¬ 
plitudes, respectively. In a superconductor, the particle- 
hole symmetry requires the zero energy spectral function 
to satisfy 39 


CA{E = O)^- 1 = A(E = 0), (1) 

where C = r y a y K is the particle-hole operator, with K 
the complex conjugate operator and r Xi y iZ the Pauli ma¬ 
trices in particle-hole space. A straightforward calcula- 
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tion shows that 


off /V— 1 


CA on C 


^ A oSy. 

0 (- d 0 a 0 + d- a) 


1 y u y\ 


lyUy 


[(-d 0 a 0 + d • a)] 1 


( 2 ) 


Comparing Eq. ( 2 ) with Eq. ( 1 ), we find that do must 
vanish at the Fermi surface. This implies that the 
particle-hole symmetry completely forbids spin-singlet 
pairing but allows spin-triplet pairing in the zero-energy 
BOSF. For the D class TSC, when there is one MF 
71 = c K + c^, on the boundary, the Majorana condition 
7 i = 7 i requires n = ft', which results in a spin-triplet 
BOSF, A° ff = ( T X <7 X — r y (jy)/2 , by choosing k =t- If 
there is an additional zero mode, 72 = i(c K — cj.), the 
anomalous spectral function is zero because the two MFs 
can emerge as a single electron or hole 39 . For a Dill 
class TSC, although there is always at least a pair of 
MFs on the boundary, the time-reversal symmetry guar¬ 
antees that these two MFs belong to different spin chan¬ 
nels and therefore manifest spin-triplet pairing 39 . The 
above analysis applies equally well to the multi-MF situ¬ 
ation and we conclude that an odd number (of pairs) of 
MFs in the D (Dill) class have only spin-triplet super¬ 
conducting correlations. It is noted that our symmetry 
analysis is independent of the details of TSC models and 
robust against the spin-independent disorder. Therefore 
the spin-triplet superconducting correlation is a universal 
and robust property of MFs. 

MF-induced resonant spin-triplet Andreev reflec¬ 
tions The current experimental methods for detecting 
MFs 17-20,22 are related to ARs. Therefore, to experimen¬ 
tally confirm our prediction of the universal spin-triplet 
superconducting correlations of MFs, we first focus on 
the AR in spin-singlet and spin-triplet channels which 
we discuss below. The scattering matrix at the interface 
of a one-dimensional (ID) SC/NM junction generally has 
the form 40 ’ 82 


R(E) = 


r ee (E) r eh (E ) 

fhe(F0 nih(£0 


r eh (E) = a 0 cr 0 + aj(jj = -a y rl e (-E)a y , 
r e e{E) = b 0 (To + bjdj = a y rl h (-E)(jy, 


(3) 


where r ( ee? hh) and ^eh,he are the normal reflection and AR 
matrices respectively, and 6 o,x,y,^ and ao, x ,y,z are the nor¬ 
mal reflection and AR coefficients respectively. The AR 
in the channel ao always couples the electron and hole 
with opposite spin, which is independent of the direction 
of the spin quantization axis 39 , and therefore provides the 
spin-singlet channel. On the other hand, we can always 
find a spin basis, where the AR in one of the channels 
a x,y,z couples the electron and hole with the same spin 39 
and therefore belongs to spin-triplet channels. As MFs 
have only spin-triplet superconducting correlations, we 
expect that the MF-induced resonant AR can occur only 
in spin-triplet channels. Below, we first use the scattering 
matrix theory 40,S2 to show the validity of our expectation 
at E = 0 and then numerically explore the AR for the 



FIG. 1. (Color online) (a) Time-reversal broken TSC/NM 
junction. The black arrows represent the spin-singlet bulk 
superconducting gap. The green arrows represent the MF 
superconducting correlations. The green curves indicate the 
spatial distribution of MFs. The red peak indicates the lo¬ 
cation of the barrier with the potential U = 2t. (b) (c) The 
superconductor is in the topological nontrivial regime. The 
conductance and AR probabilities in spin-singlet and spin- 
triplet channels as a function of the energy inside the su¬ 
perconducting gap. (d), (e) The superconductor is in the 
topological trivial regime. The conductance and AR proba¬ 
bilities in spin-singlet and spin-triplet channels as a function 
of the energy inside the superconducting gap. In (c) and (e), 
the solid blue, solid red, dashed black and dashed green lines 
represent AR probabilities, |ao| 2 , \a x \ 2 , \a y \ 2 , and \a z \ 2 for 
spin-singlet and spin-triplet channels respectively. 


energy range inside the superconducting gap. For the D 
class TSC/NM interface, it was pointed out 42 that the 
incoming electron and reflected hole in the MF-induced 
resonant AR process have the same spin. Therefore, this 
AR can occur only in spin-triplet channels. For a Dill 
class TSC/NM interface, the condition TR^T -1 = R, 
due to time-reversal symmetry, constrains b X:V:Z = 0. 
Here, T = ia y K. The unitarity of R further requires 
either (1 60 1 2 + |ao| 2 = 1 and a x = a y = a z = 0) or 
(\a x \ 2 + \a y \ 2 + \a z \ 2 = 1 and b 0 = a 0 = 0) 39,S1 . The 
former corresponds to the trivial spin-singlet AR with 

Q = Pf ^ i(j y R ^ = l 40 ’ S2 . More importantly, the lat¬ 
ter indicates a perfect spin-triplet AR and results in the 
nontrivial topological invariant Q = — l 40 ’ S2 . 

To further explore the AR inside the entire super¬ 
conducting gap, we perform numerical calculations of 
ARs in two concrete ID TSC tight-binding models. The 
Hamiltonian of the first TSC model, which violates time- 
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reversal symmetry, has the form 4,5 

Hts = [—2t cos(fc) - fi s ] r z (8 cr 0 Mr 0 (8) cr z 

+ 2t so sin(fc)r 2 (8 cr y + A r x (8 cr 0 , (4) 

where t is the spin independent hopping, t so is the SOC 
hopping and fi s is the chemical potential of the SOC wire 
with M the Zeeman coupling strength and A the prox¬ 
imity induced superconducting gap. The Hamiltonian of 
the normal metal is 

-Hnm = (-2 1 cos(fc) - n n )r z <g> a 0 , (5) 

where /i n is the chemical potential. We numerically calcu¬ 
late the scattering matrix at x = x\ (shown in Fig. 1(a)) 
based on the Fisher-Lee relation 44 , 

Rij(E,r) = -Sij + ihy/ViVjGfj{E,r), (6) 

where G R is the retarded Green’s function with = 
1,..., 4 and is the velocity of the particle at energy 
E in i(j) channels. In Fig. 1(b), by adding a barrier 
potential U = 2t at the TSC/NM interface and taking 
M = 0.8£, A = t s o = 0.1 £, we plot the calculated con¬ 
ductance as a function of energy. The quantized zero- 
bias conductance with height 2 e 2 /h confirms the exis¬ 
tence of a MF at the TSC/NM interface. Moreover, we 
plot the absolute value of the AR probabilities \ao iXj y iZ \ 2 
in Fig. 1(c). It is shown that inside the gap, the AR only 
occurs in the spin-triplet channels a x and a y , which is 
consistent with our prediction. In Figs. 1(d) and 1(e), we 
also plot the conductance and the AR probabilities when 
the superconductor is in the trivial regime by choosing 
M = A/2 = 0.05L The conductance has no peak at 
E = 0 and the AR only occurs in the spin-singlet chan¬ 
nel. These results suggest that the perfect AR in the 
spin-triplet channels is solely induced by the MF. 

Our prediction of the spin-triplet MF superconducting 
correlations is universal for both time-reversal violated 
and respected superconductors. Therefore, in the fol¬ 
lowing, we choose a time-reversal symmetric TSC model 
whose Hamiltonian has the form 9 

Hts = [—2 t cos(fc) - /i s ] r z (8 ao 

+2t so sin(fc)r z (8 o z + A(k)r x (8 cr 0 , (7) 

where A(fc) = [Ao — Ai cos(fc)], and is closed at 
cos(fco) = Aq/Ai. The bulk superconducting gap in 
Eq. (7) is a combined s± pairing potential and is ob¬ 
viously spin-singlet. The SOC in the semiconductor 
nanowire induces the spin splitting which results in two 
Fermi wave vectors fcif and & 2 f. We assume fcif < fc 2 f, 
as shown in the inset of Fig. 2(a). It was pointed out 9 
that for fcif < fco < fc 2 f, the system is in the topologi¬ 
cal superconducting regime. In the following, we choose 
Ai = 2A 0 = 0.2£, fi s = —£, and t so = 0.3t. The gap 
sign is negative at k = fcif and positive at fc = fc 2 f as 
shown in Fig. 2(a). We numerically calculate the reflec¬ 
tion matrix at x\ (Fig. 2(a)). The conductance extracted 
from the reflection matrix is plotted as a function of en¬ 
ergy in Fig. 2(b). The quantized 4e 2 /fc zero bias-peak 



FIG. 2. (Color online) (a) Time-reversal invariant TSC/NM 
junction. The white arrows represent the spin-triplet super¬ 
conducting correlations of MFs. The black arrows represent 
the spin-singlet bulk superconducting gap. The bulk disper¬ 
sion of the TSC shows two superconducting gaps. The inset 
shows the sign of the minimal gap A m i n is negative, (b) The 
conductance of the TSC/NM junction as a function of energy, 
(c) AR probabilities for spin-singlet and spin-triple channels 
as a function energy. The solid blue, solid red, dashed green 
and dashed black lines represent AR probabilities, |ao| 2 , \a x \ 2 , 
\a y \ 2 and \a z \ 2 respectively. 


indicates that there is a pair of MFs at the end of the 
time-reversal invariant TSC. More importantly, we plot 
the AR probabilities lao^,^! 2 a t x i as a function of the 
energy in Fig. 2(c). Inside the gap, we find that the AR 
occurs only in the spin-triplet channel a z . The underly¬ 
ing physics is that the superconducting gaps at Tfcff and 
±fc 2 f have opposite signs. Therefore the AR coefficient 
associated with Tfcif also has the opposite sign to that 
associated with ±fc 2 f and thus the AR matrix is in the 
spin-triplet channel with the form a z a z instead of the 
spin-singlet channel with the form aocro- Therefore, al¬ 
though the bulk superconducting gap is spin-singlet, the 
AR at the TSC/NM interface is in the pure spin-triplet 
channel. 

SOC controlled critical Josephson current- So far, we 
have shown that the MF-induced resonant AR occurs 
only in spin-triplet channels. It is known that the spin 
property of the AR-induced Cooper pairs is the same 
with that of the associated AR 45 . Therefore, the MF- 
induced resonant AR always injects spin-triplet Cooper 
pairs. Recent studies 45-48 show that SOC can rotate the 
d vector of the AR-induced spin-triplet Cooper pairs but 
leave the spin-singlet pairs unaffected in the proximity re¬ 
gion. Thus SOC provides a mechanism to distinguish be¬ 
tween spin-triplet and spin-singlet Cooper pairs. To con¬ 
firm our prediction, we therefore consider a TSC/SOC- 
semiconductor/TSC Josephson junction (Fig. 3(a)). The 
Hamiltonian of the SOC-semiconductor system has the 
form 

#semi = [-2i cos(fc) - fj, B ] t z ® <r 0 + 2i' 0 sin (k)r z <g> a y , 
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FIG. 3. (Color online) (a) TSC/SOC/TSC junction, (b) An¬ 
dreev levels for different SOC strengths which correspond to 
the precession angles 0 = 0, 37r/4, 7 t, 57t/4, 2tt, respectively, 
(c) The oscillatory critical current as a function of 0. The 
current is normalized by the critical current at t' so = 0. 


where t' so is the SOC hopping amplitude and the effective 
magnetic field of the SOC is along the y direction. On 
the boundary x\ and X 2 [Fig. 2(a)] of time-reversal bro¬ 
ken TSCs, based on our previous discussion, the Cooper 
pairs can be described as | tt)> with spin t along the 
z direction. When a Cooper pair travels from x\ to X 2 
in the SOC-semiconductor region, because the effective 
SOC field is along the y direction, the spins of the elec¬ 
tron and hole start to precess in the x-z plane with the 
same precession angle. As the spin splitting energy due 
to SOC is SE = 2t' 0 sin(fcf), the precession angle of the 
particle traveling from x\ to X 2 satisfies 39 


<9 = 


5EL 
hv f 


fj 

= TV—, 
t 


( 8 ) 


where Vf ~ 2£sin(/q) in the limit t' so £, and we assume 
there are N sites in the normal SOC region whose length 
is L (Fig. 3(a)). The state of this Cooper pair at X 2 
generally can be described as | /*) where represents 

a spin lying in the x — z plane and its spin wave function 
is (cos(0/2), sin(#/2)) T . Therefore, the projection of the 
condensate | f) to | tt) takes the form 


(tt I /*) = (l + cos#)/2. (9) 


Especially for 6 = 7r, the projection in Eq. (9) is zero 
which means no Cooper pair can travel through the 
Josephson junction, and this indicates zero critical cur¬ 
rent. As a result, we expect the critical current will oscil¬ 
late with a large amplitude and may even vanish as the 
SOC strength is suitably tuned. To confirm our theoret¬ 
ical argument, we plot the Andreev levels in Fig. 3(b) in 
the cases of 0 = 0,37r/4,7r, 57t/4, 2tt respectively. Because 
the Andreev levels cross E = 0 at </> = 7r, the Josephson 
current has the form I s = I c sin(0/2) where I c is the crit¬ 
ical current and </> is the phase difference between two 
TSCs. It is noted that at 0 = 7r, the energy-phase rela¬ 
tion of the Andreev levels is almost fiat. Therefore we 
conclude I c = 0 in this case which is consistent with 
our theoretical analysis. In Fig. 3(c), I c is plotted as a 
function of SOC controlled precession angle #, showing 
a 27r period as we expected based on Eq. 9. Generally, 
time-reversal broken TSCs can be considered as Cooper 
pair filters which only allow the Cooper pairs in the state 
| tt) to enter the TSCs from the normal region. However 
the situation for time-reversal invariant TSCs is very dif¬ 
ferent. Because the superconducting condensates exist 
in both spin-up and spin-down channels, the spin-triplet 
Cooper pairs can always enter the time-reversal invariant 
TSCs which is independent of their spin direction 39 . 

Our prediction can also be applied to multi-band 
and even two-dimensional quantum wells with equal 
Rashba and Dresselhaus SOC strengths, which have 
been realized experimentally 49-54 . In this case, the 
spin precession angles for all transverse channels are 
the same 55-57 and the d-vector of the triple pairs will 
not decay even in the presence of spin-independent 
disorder 45 . The oscillatory critical current has also been 
predicted and observed in the non-topological super- 
conductor/ferromagnet/superconductor (SFS) Joseph¬ 
son junction where time-reversal is broken and MFs do 
not play any role in contrast to our system 58,59 . The 
current-phase relation in the SFS junction undergoes 
a 0-7T transition when the critical current approaches 
zero. We emphasize that the current-phase relation of 
our time-reversal broken Josephson junction does not un¬ 
dergo a 0-7T transition. This makes our prediction quali¬ 
tatively different from that in the SFS junction. There¬ 
fore, our prediction of a SOC controlled critical current 
is a transport signature of MFs. 

In a realistic experimental setup, the semiconductor 
wire may not necessarily be perfect ID but may have 
several channels. However, our prediction requires the 
scattering at the TSC/NM interface to be ID. The re¬ 
cent development of a surface split-gates design 60 can 
make a real ID point contact between the topological su¬ 
perconducting semiconductor wire and the normal lead. 
Therefore, our experimental proposal can be realized by 
an existing experimental technique. 
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Supplementary information 


Appendix A: Some knowledges 
1. Particle-hole symmetry and time-reversal symmetry 

The BdG Hamiltonian, H , of a superconductor has always particle-hole symmetry CHC = — H where C is the 
particle-hole operator. In the basis (c^, , c j, —cj) T , the particle-hole operator is C = r y a y K where r and a are Pauli 

matrices in particle-hole and spin space respectively. The time-reversal operator has the form T = irocr y K. If the 
system is time-reversal invariant, we have T~ 1 HT = H. 


2. Definition of spin-singlet and spin-triplet gaps 


Theoretically speaking, the Hamiltonian for superconducting pairing potential generally has the form 

Hgap = (AoTo + Aofo) + (A Ayfj), 

i= 1 , 2,3 


(Al) 


where Ao(Aq) is the real (imaginary) part of the s-wave gap, Aj(Aj) is the imaginary (real) part of the p-wave gap, 
To,1,2,3 = ( T x ®a 0 ,Ty®(T x ,Ty®<Jv,Ty®<Jz ) and f 0 ,i,2,3 = ( Ty®(jQ,Tx®(T x ,T x ®<jy,T x ®<j z ) in the basis (c t , cj., cj,-c|) T . 

We find that the matrices To, 1,2,3 and To, 1,2,3 satisfy the Dirac algebra 

{r M ,W} = 2/4x4, jf M ,fV j = 2 74x4, f^o,1,2,3 = ^5^=0,1,2,3, (A2) 

where T 5 = irir2r 3 r4 = —ir z ao. Therefore, the pairing potential can be constructed from these eight gamma 
matrices. In the following, we will prove that T*-4,2,3 expand a three dimensional vector space and construct a three 
dimensional representation of a SU(2) group whose generators correspond to the rotation in spin space. 

First, we define 

= '^ijk T 0 ® G k — 2 e ijk ^ q ^ = e ijkSki (A3) 


where £*=1,2,3 are three generators of a SU( 2 ) algebra. Next, we will show that the matrix e lS ri corresponds to 
rotating the spin from z axis to the n direction where n is a three dimensional unit vector. Generally a wave function 
in the basis (c^, c^, cj, — c|) T has the form 


T = 


iCTy'lph 


(A 4 ) 


where ^ e (h) is the 2 -component spinor wave function for the electron (hole). When we rotate the electron spin as 

V/ = e— / 2 V>e, (AS) 

the spinor wave function for the hole will be transformed as 

< = e <<r *-"/ Vh. (A 6 ) 


Therefore the wave function is transformed as 


( € \ = ( e-™' n rp e 

V Wy-K ) V 


2 ip e 

e~ ia " n ia y ip h 


— e ~ iS n ^>. 


(A 7 ) 


Therefore the matrix e lS n corresponds to the rotation in the spin space. Because the SU(2) generators Si— 1,2,3 
are constructed from the Eq. A 3 , the three gamma matrices Ti, 2,3 construct a three dimensional representation of 
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a SU(2). Therefore Ti 5 2,3 are the spinnors of the SU(2) with generators Sk = This is consistent to the 

fact the T 12,3 are the basis of spin-triplet pairing. It is easy to see that (Ti,^) = {ir y a x ,ir y a y ) represent the 
superconducting pairing in the same spin channel. Therefore, as long as the superconducting pairing can be written 
as the superposition of the basis Ti ? 2 , 3 , we can always find a rotation operator which can rotate a general spinor to 
the axis in which the electron and hole have the spin. Also it is easy to check that [T 0 , ^= 1 , 2 , 3 ] = 0. Therefore T 0 is 
a one-dimensional representation of the SU(2) spin rotation group which is consistent to the fact that T 0 is the basis 
of the spin-singlet pairing. The similar results can be applied to T 0 , 1 , 2 , 3 - 


Appendix B: Spectral function and pairing density 


The equal space spectral function at energy E is defined as 

A(E,x) = ±-(G r (E,x) - G a (E,x )), 

which is a Hermitian matrix. Based on the Lehmann representation, the retarded Green’s function is related to the 
spectral function as 

G R {E,x) = (E-H + iS)~ 1 = f _EE> de = P f EE. de - inA(E,x), 

J -oo E -e + id J _ 00 E - e 

where x and x' are the spatial coordinates. Because the spectral function A(e, x) is hermitian, the principle 
integral P is hermitian, and —inA(E, x) is anti-hermitian. The spectral function can be written as 

^(e) = Ee„=J<M x))(T n (x)|. As we are interested in the superconducting pairing, in the following, we will fo¬ 
cus on the off diagonal part of the spectral function which has the form A off (e) = d /x (e,x)T M + d M (e, x)Y^ where d y 
and d y are the superconducting pairing amplitudes. Because T M , T M and the spectral function A(e,x) are hermitian, 
the coefficients d M and d y are real. We now focus on the retarded Green’s function at E = 0 because it is related to 
the MFs in the topological superconductors. As the Hamiltonian H has particle-hole symmetry, the retarded Green’s 
function at E = 0 also respects particle-hole symmetry which is shown below: 

CG r (0,x)C~ 1 =C(-H+ iS)- 1 C~ 1 =-{-H+ IS)- 1 =-P f Eh. x b e + inA(E, x) (Bl) 

J — OO ^ 

Particularly we are interested in the block off diagonal part of the retarded Green’s function 

fR =(7“ C) - <B2) 


Based on Eq. Bl, the off diagonal retarded Green’s function satisfies 

CF r ( 0, x)C~ 1 = P f EEE de + *7rA off (0, x) 

J — 00 e 


= p 


/ 


T d^(e, x)Fp 


(d^0,x)T^+d^0,x)f, 


On the other hand, particle-hole symmetry requires the anomalous Green’s function to satisfy 


CF K (0,x)C- L = C P 


•/: 


A off (e,a;) 


—e 


— i7rA off (0, x) J C~ x 


= -P 


/: 


d 0 (e, aj)r 0 + d 0 (e, :r)f 0 - d,(e, x)Yi - di(e,x)ti 


(B3) 


+ 


in do(0, x)r 0 - d 0 { 0, x)t 0 + d t ( 0, x)T i: + ^(0, x)t^j , 


(B4) 


where we use the fact that 

CT M= o,1,2,3^ = ( T x (jQ,T y (7 x ,T y (7 y ,T y (J z ), CY yi= 0,1,2,3^ = ( Ty<To, T X (7 X , T x (7 y , T X (7 Z ). 

Comparing Eq. B3 with Eq. B4, we conclude that spin-singlet pairing is completely forbidden in the anti-hermitian 
part of the anomalous Green’s function F^ rlti _ herrnitian (0, x) = —i7rA off (0,x) which indicates that the MFs cannot 
have spin-singlet pairing. 



Appendix C: Scattering matrix theory for Dill topological superconductor 


When a normal lead is attached to a superconductor, the Andreev reflection at the SC/NM interface will induce the 
superconducting condensate in the normal region. In the following and for a one-dimensional (ID) SC/NM junction, 
we use scattering matrix theory to show that the induced superconducting condensate by the MF resonant Andreev 
reflection is always spin-triplet. We first focus on a Dill class SC/NM junction. Because of time-reversal symmetry, 
the scattering matrix is self dual si which requires that 


S = 


T'ee&O CL0&0 T 

aa 0 - djcjj r hh ao 


It is straightforward to show that the above scattering matrix satisfies the condition T (E)f- 1 = S(E). It is also 
noted that aj =X:VrZ are the spin-triplet scattering amplitudes and ao is the spin-singlet scattering amplitude. We are 
interested in the Andreev reflection at E = 0. In this case, according to the particle-hole symmetry of the system, we 
have 


7~y(Jy 


= s = 


0 -i<jy \ / r ee cro ao&o + ajCTj 

icjy 0 ) \ aa 0 - djcjj r hh a 0 

^hh^O ~( a 0 cr 0 + a pj) 

-(a*cr 0 - a*CFj) r* e <7o 

Tee^o a^CTo CLjCTj 
a 0 (T0 QjjCJj ^hh^O 


0 — i(Jy \ 

i(Jy 0 J 


Therefore, we have r* e = r^h = fro is real, ao and aj are pure imaginary. At last, the matrix S should be unitary. 
This requires 


S^S = 



0 


Because 


S) ee — (|fro| 2 + |^o| 2 + \d x \ 2 + \(iy\ 2 + l^l 2 )^ + 2aoaj<jj — ao, 

we have either aj =x ^ y ^ z = 0, 16o| 2 + |ao| 2 = 1 or ao = 0, 16o| 2 + \a x \ 2 + \a y \ 2 + \a z \ 2 = 1. These indicate that the 
spin-singlet Andreev reflection (the former case) and spin-triplet Andreev reflection (the latter case) cannot coexist. 
For the spin-triplet Andreev reflection, we further have 

(S^S) e h = 2 boaj<Tj = 0, 

which requires either r* = 0 or aj =X:VrZ = 0 and indicates that the normal reflection and spin-triplet Andreev reflection 
also cannot coexist. Therefore, we conclude for a Dill class SC/NM junction and at E = 0, the spin-triplet Andreev 
reflection must be a perfect Andreev reflection with vanished spin-singlet component. It is also straightforward to 
show that 


Pf ia y r 0 


0 

~ a j (T j 


= Pf 


0 i ct yCtj a j 


% CT y CL j CT j 

|2 i 12 


0 


— —{\d X \ + \dy\ + \a z \ ) — —1, 


which is consistent to the condition of the nontrivial topological quantum number S2 in terms of the scattering matrix 
for a Dill class superconductor. 


Appendix D: Spin precession angle and its effect for time-reversal invariant topological junction 

The Hamiltonian of the SOC-semiconductor is 

Hsemi = (—2tcos(fca) - fi s )r z 0 a 0 + 2t' 0 sin (ka)r z 0 a y , (Dl) 

where k is the wave vector, a is the lattice constant, t' so is the SOC hopping constant in the SOC-semiconductor and 
the effective magnetic field of the SOC is along y direction. SOC will shift the two spin-subbands to opposite direction 
in the momentum space. In the ID tight-binding model with k x a y -\ike SOC Hamiltonian at low energy limit, the 
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dispersion relation is E = 2t(l — cos (ka)) ± 2t so sin (ka) = 2t — 2^/t 2 + t 2 0 cos(fca db 50) where cos(50) = t/^/t 2 + £g 0 • ^ 
there are N sites in the normal region, the electron and hole will take time T = N/vf ~ Nh/2at sin(&/a) to go through 
the normal region. Because the effective magnetic field due to SOC is B soc = 2t' 0 sin(fc/a), the spin precession angle 
of an electron or hole in the normal region is 7 = B soc T = Nt' so /t. In the continuous limit, we have t —> h 2 /2ma 2 , 
t ' s o a/a. Therefore, the spin precession angle in the continuous limit is 7 = Nt' so /t = 2maL where L = Na is the 
length of the normal region. 

For a time-reversal invariant TSC, there is a pair of Majorana zero modes with opposite spins at the boundary. 
The induced spin-triplet condensates are described by (| tt) + | 00))/\/2 which is even under time-reversal operation. 
Here we assume the Majorana spins are parallel or anti-parallel to z axis. If the spin-orbit coupling (SOC) field 
direction is along y axis as we consider in the manuscript, the spins are rotated so that the spin-triplet states are 
changed to (| /*/*) + | v / v /))/\/2 where and represent the spins lying in the x-z plane with wave functions 
(cos(0/2), sin(#)/2) T and (sin(0/2), — cos(#)/2) T respectively. Therefore the projection of the condensates is 

((tti + aii)(i^) + i/^))/2 = i, 

which is independent of the SOC. Therefore in this case the TSC/SOC/TSC junction performs like a normal Josephson 
junction. However if the SOC field direction is along x axis so that the Majorana spin is rotated in y-z direction, the 
spins and correspond to (cos(0/2), i sin(#)/2) T and (i sin(0/2), cos(6 > )/2) T respectively. The projection is 

((tt I + (ii |)(| /*/*) + \)/ v/))/ 2 = cos0, 
which will result in a 0- and tv- Josephson junction transition. 
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